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I. INTRODUCTION 



G. Pagliara, J. Schaffner-Bielich 
Institut fiir Theoretische Physik, Goethe Universitdt, D-60438, Frankfurt am Main, Germany 

We study the equation state of strongly interacting quark matter within a NJL-like model in 
which the chiral condensates and the color superconducting gaps are computed self-consistently as 
a function of the baryon density. A vector interaction term is added to the Lagrangian in order 
to render the quark matter equation of state stiffer. For the low density hadronic phase we use 
a relativistic mean field model. The phase transition to quark matter is computed by a Maxwell 
construction. We show that stable CFL cores in hybrid stars are possible if the superconducting gap 
is sufficiently large. Moreover we find stable stellar configurations in which two phase transitions 
occur, a first transition from hadronic matter to 2SC quark matter and a second transition from 
I 2SC quark matter to CFL quark matter. 
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The possibility that quark matter, and eventually color superconducting quark matter, is present in the center of 
neutron stars has stimulated many theoretical investigations in the last years both on the modelling of the equation of 
state (EoS) of quark matter and on the phenomenological signatures of the presence of quark matter in neutron stars 
(~| I P] . Presently the "state of the art" for the EoS of quark matter is represented by the three-flavor Nambu- Jona-Lasinio 
Qh' (NJL) model in which both the chiral condensates and the diquark condensates are self-consistently computed as a 
function of the chemical potential and temperature. Chemical equilibrium and charge neutrality (both electric and 
color charge) conditions necessary to describe neutron star matter are also imposed in this model [1, H, 0, S @]- 
^ ' Recently, also the effect of a finite neutrino chemical potential has been included 0, @|- The structure of the QCD 
d . phase diagram within this model turns out to be very rich, with many different possible quark phases. One of the 
most striking feature, on which we will focus here, is the first order phase transition between the two flavor phase, 2SC 
or normal quark matter depending on the diquark coupling constant, and the three-flavor superconducting phase, the 
J> ' Color-Flavor-Locking phase (CFL), at vanishing temperature. 

[ Concerning the phenomenological signatures, the differences between the mass-radius relation for neutron stars 
■ and quark or hybrid stars is currently studied. The astrophysical data so far are still affected by large uncertainties 
but will improve considerably with the advent of new satellite missions, as XEUS, Constellation-X, SKA, JWST and 
" I . LISA. From recent theoretical studies it turns out that the values of the maximum mass of neutron stars and compact 
^-H ' stars containing quark matter are very similar [lOl. Therefore, it seems difflcult with the present knowledge to 
rule out quark matter from just a mass measurement [ill. [T^. Other interesting quantities have been calculated (see 
Ref- [IS] for a recent review) for the different possible phases of quark matter, as the neutrino emissivity and the 
heat capacity which are important for the cooling of compact stars [13, [iM H^ - ij . 1 ^ or the bulk viscosity which 



K>" ■ determines stability with respect to gravitational waves emission via r- modes [Tgf. l20l [2iri22l. [23l [2^ . Also in explosive 
Kyi ' phenomena, as supernovae and gamma-ray-bursts (GRBs), quark matter can play an important role. For instance, 
. the possibility of a double phase transition, first from hadronic matter to 2-flavor quark matter and then from 2-flavor 
' quark matter to the CFL phase has be en prop osed to explain the complicated time structure of GRBs exhibiting a 
long quiescent time in their light curve [25l [261 . [27| . 

When the above mentioned NJL- EoS is used for the applications on compact stars, hybrid stars become unstable at 
the onset of the CFL phase and therefore CFL phase can not be present in the core of neutron stars. This conclusion 
was obtained first in Ref. [1^ where quark matter does not occur at all in compact stars because there is a direct 
transition from hadronic matter to CFL matter in the model used. In Refs. [l^, [13, [Ml for different EoSs for the 
hadronic matter and different parameters for the NJL model for quark matter, again the CFL phase was ruled out 
because it renders the stars unstable. The 2SC phase could appear; the conclusion is therefore that only 2-flavor 
superconducting quark matter can be realized in compact stars. A similar result about the allowed quark phases in 
compact stars was also found in Ref. j32] where a modified 2-fiavor NJL model, which simulates confinement at low 
density, is proposed. The dependence of the stability of a quark core from the momentum cut-off of the NJL model 
has been analysed using density dependent cut-offs, but the instability still persists [s^ [6l| . 

On the other hand, completely different results are obtained using MIT-bag-like models as shown in Refs. [3^. [ssl. [36j 
where the appearance of CFL cores does not compromise the stability of the star. Moreover, the absence of the 2SC 
phase in compact stars has been demonstrated in Ref. . 

In this paper we want to consider again the NJL- EoS as computed in Ref. [1] to study the structure and composition 
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of compact stars. We will investigate larger windows of the model parameters with respect to previous work, with 
particular attention to the diquark coupling. Furthermore, we discuss the importance of the procedure used to fix 
the effective bag constant within the NJL model for the stability of a star when the phase transition to quark matter 
is considered. We will investigate a new procedure for fixing the effective bag constant by requiring that the chiral 
symmetry restoration coincides with the transition from the hadronic to the quark matter description [6^ . We use 
then different EoSs to compute the mass-radius relations of hybrid stars showing that in some cases a stable CFL 
core is possible and, even more intriguing, that two phase transitions, from hadronic matter to the 2SC and then to 
the CFL phase, can take place in compact stars. The double phase transition is particular intriguing in connection 
with the deconfinement quark model of GRBs in which the interpretation of bursts presenting two emission periods 
is due to a double phase transition in compact stars [H, [2^, . 

The paper is organised as follows. In Sec. II we discuss a simple toy model for a first order phase transition between 
hadronic matter and quark matter in compact stars and we study the stability of the star by varying the parameters 
of the quark matter EoS. In Sec. Ill we compute the quark matter EoS within the NJL model for different sets of 
parameters and finally in Sec. IV we discuss the stability of the stars obtained using our EoSs with particular emphasis 
on the CFL core stability. In Sec. V we draw our conclusions. 

II. A TOY-MODEL FOR PHASE TRANSITIONS IN COMPACT STARS 

We present in this section a toy-model EoS for strongly interacting matter to show qualitatively the conditions for 
stable solutions of the Tolman-Oppenheimer-Volkoff (TOV) equation when a phase transition from hadronic matter to 
an exotic phase, like quark matter, occurs. At very low baryon densities, n < O.OSfm^'^, we use the EoS as computed 
in Ref. [38[ (similar results are obtained using the EoS of Ref. [s^) suitable for the crust of neutron stars. At larger 
densities, we consider for the hadronic matter a relativistic mean field EoS, GM3, taken from Ref. [i^. For quark 
matter we adopt a schematic EoS in which the pressure is proportional to the energy density p = ae where the slope 
(which corresponds to the sound velocity and therefore regulates the stiffness of the quark matter EoS) is a free 
parameter. Notice that this EoS corresponds, for a = 1/3, to the EoS of massless and non-interacting quarks. We 
then model the phase transition from hadronic matter to quark matter by introducing two other free parameters, 
the energy density jump = 62 — ei at the onset of the phase transition (e2 and ei are the energy densities of the 
quark phase and of the hadronic phase) and the pressure at which the phase transition occurs po- Our aim is to 
investigate, for the different values of these parameters, the stability of stars having a core of quark matter. Let PI be 
the parameter set with a = 1/3, Ag = 200 MeV/fm"^ andpo = 200 MeV/fm^. In panel (a) of Fig.[T]the corresponding 
EoS is plotted by a gray line. In this case compact stars are unstable when the quark matter phase appears, see the 
gray line in panel (b). If we decrease now the transition pressure to po — 50 MeV/fm'^ keeping the same values for 
the other parameters (set P2 see solid line in panel (a)), a sizable branch of stable configurations appears instead, 
see the thin solid line in panel (b). If we reduce the energy density jump to A^ = 50 MeV/fm'^ (set P3 with po and 
a as in set PI, see dashed line in panel (c)) a small branch of stable solutions is obtained (dashed line in panel (d)). 
Concerning the stiffness parameter a, to obtain stable solutions one must choose a = 1 (set P4 with po and A^ as 
in PI). For this value of a the quark matter EoS is even stiffer than the hadronic matter EoS (see the dotted line 
in panel (e) for the EoS and the dotted line in panel (f) for the mass-radius relation). This last case is probably not 
realistic in view of the asymptotic freedom property of QCD. 

In conclusion, what seems to be the most crucial parameter for the stability of a quark matter core is the value of 
the pressure at the onset of the phase transition. If the softening due to the appearance of quark matter occurs at too 
large pressures, deep in the core of a neutron star, the sourronding material exerts a pressure that can not be sustained 
by the new formed phase and therefore the star collapses. Also a low value of A^ can help to stabilise the quark 
matter core. A detailed analytical study on the critical value of energy density jump for having stable cores of a new 
phase can be found in Refs. [4l|, li^, IH, Q demonstrating that £2/61 should be not larger than £2/61 < |(1 +Po/^i)- 

We will see in the next sections how these cases are connected with more physical meaningful quantities as the bag 
constant, the superconducting gap and the constituent masses of quarks . 

III. PHASE TRANSITION TO QUARK MATTER 

We present now a more realistic EoS for quark matter. The EoS is computed within the NJL-like model proposed 
in Ref. :3"] in which a scalar diquark interaction term for the color antitriplet and flavor antitriplet channel is added 
to the usual NJL model. Here we include also the isoscalar vector term as in Ref. [13, IH, Hg] in order to obtain 
stiffer EoSs. The input variables of the model are the chemical potentials for all the quark flavors and colors given, 
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FIG. 1: Equations of state and corresponding mass-radius relations for the toy- model for different parameter sets. EoSs are 
shown in the left row and mass-radius relations in the right row. The purely hadronic matter EoS and the corresponding 
neutron stars sequence are depicted by a thick solid line. Panel (a): quark matter EoS for a = 1/3, Ae = 200 MeV/fm"^ and 
po = 200 MeV/fm^ (gray line set PI), and a = 1/3, = 200 MeV/fm^ and po = 50 MeV/fm^ (solid Una set P2). Panel (c): 
quark matter EoS for the choice a = 1/3, A^ = 50 MeV/fm^ and po = 200 MeV/fm^ (dashed line set P3), the other curves are 
as in panel (a). Panel (e): quark matter EoS for a = 1, A^ = 200 McV/fm^ and po = 200 McV/fm^ (dotted line set P4), the 
other curves as in panel (a). Stable configurations with a quark core can only appear for sufficiently stiff quark matter EoS or 
for a transition point away from the maximum mass configuration of the purely hadronic compact stars. 

in chemical equilibrium, by the matrix: 

= {^^s-^ + ^^QQf) Sa6 + [m3 (Ts),, + MS {Ts)j s-f" . (1) 

where jj, is the quark chemical potential, hq is the chemical potential of the electric charge equal to minus the electron 
chemical potential /ig, Ms and /itg are the color chemical potentials associated with the two mutually commuting color 
charges of the SU{3)c gauge group. The cxpHcit form of the electric charge matrix is Qf = diagj(|, — |, — 5), and for 
the color charge matrices T3 = diag^{^,—^, 0) and V^Ts = diagg(^, 5, — !)• 

The model parameters are fixed by fitting low energy hadronic properties which are the current quark masses, the 
quark- antiquark coupling Gs, the strength K of the '"t Hooft" interaction and the cut-off parameter A introduced in 
the NJL model to regularise the ultraviolet divergences: 

m„,d = 5.5 MeV, (2) 
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FIG. 2: Upper left panel: pressure as a function of the chemical potential for hadronic matter (thick line), the 2SC phase 
(dashed line) and the CFL phase (solid line) . Upper right panel: pressure as a function of the baryon density. Parameters are: 
Gd = Gs, B = Bo, Gv = 0. A direct transition from hadronic matter to the CFL phase occurs with a large baryon density 
jump ~ 0.45 fm~'^ at the onset of the phase transition. Lower left panel: pressure as a function of the chemical potential. 
Lower right panel: pressure as a function of the baryon density. Parameters are: Gd = Gs, B = B,, Gv = 0. For a lower 
value of B there is first a transition from hadronic matter to 2SC quark matter and then from 2SC quark matter to CFL quark 
matter. 



nis = 140.7 MeV , (3) 

GsA^ = L835, (4) 

KA^ = 12.36, (5) 

A = 602.3 MeV . (6) 



After fixing the masses of the up and down quarks by equal values, mu,d — 5.5 MeV, the other four parameters are 
chosen to reproduce the following four observables = 135.0 MeV, rriK = 497.7 MeV, to,;' = 957.8 MeV, and 

= 92.4 MeV. This parameter set gives m,, = 514.8 MeV \4T\ . 

There are two more parameters, the diquark coupling Gd and the vector current coupling Gy which are not 
known. We will use Gd = Gs and Gd — l-2G's because one expects that the diquark coupUng has the strength as 
the quark- antiquark coupling. For Gy we choose the cases Gv = and Gv = 0.2Gs- 

At vanishing temperature and within the mean field approximation the pressure reads: 

. 18 M -,3 3 2 

P = / dfcfc^|6,|+4i^a„a,a,-^^|A,|2-2G5^a^ + ^+Pe (7) 

i—l ^ ^ c— 1 a— 1 ^ 

where ei are the dispersion relations as computed in Ref. (Ju,d,s are the quark- antiquark condensates and Ac are 
the three diquark condensates. We denote with ujq — 2Gv (Q-^l V'lV'u + V'dV'd + '4'l'4's\QM) the mean field expectation 
value of the scalar vector meson lo. This field modify also the chemical potentials: ^u,d,s fJ'u,d.s ~ ^o- Finally, the 
contribution to the pressure of electrons is pe — (YI-k^^. 

The pressure within the NJL model is defined but for a constant B, similarly to the MIT bag constant, which is 
usually fixed by the following procedure: one requires that the corrected pressure p — B is vanishing at vanishing 
chemical potential |4|, i4S,, J£] . In our model, for the parameters set used here, we have B = Bq = (425.4MeV")'* }63j . 
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Actually, this procedure to determine the bag constant is somehow unsatisfying, as also stated in Ref . [49[ , since the 
pressure computed within the NJL model at vanishing density is used i.e. in a regime where NJL model can not 
be trusted due to its lack of confinement. On the other hand, in the MIT bag model for instance, which contains 
confinement, the pressure in the vacuum is not vanishing. 

Here we propose an alternative procedure to fix the bag constant in the NJL model. First, we introduce at low 
density an EoS having hadronic degrees of freedom, like the GM3 EoS used in Sec. II, and then we compute the 
transition to quark matter, the deconfinement transition, by a Maxwell construction. We remark that " deconfinement" 
in our scheme has the meaning of a change of degrees of freedom and the corresponding Lagrangian, it is not a phase 
transition described by an order parameter. To fix the bag constant we assume that deconfinement occurs at the same 
chemical potential as the chiral phase transition i.e. when chiral symmetry is restored. Practically this means that 
we require that the pressure of quark matter, p — B, is equal to the pressure of the hadronic matter at the critical 
chemical potential for which chiral symmetry is restored i.e. the value ^icrit computed in the NJL model. This allows 
to fix the value of B and, as we will see, to obtain significantly different results with respect to the ones obtained using 
the conventional procedure. It turns out that deconfinement occurs at very large quark chemical potential ~ 470 MeV 
by far larger than the critical chemical potential for chiral symmetry restoration for B — Bq, the standard choice [49l |. 
The bag value obtained with our assumption, B — B,,, as we will see in the following is marginally smaller than Bq 
and must be considered as the lowest possible value for the bag constant in the NJL model because it allows to use the 
NJL-EoS just starting from ficrit- For chemical potentials lower than ficrit the density of quarks, as computed within 
the NJL model, is vanishing due to the completely broken chiral symmetry. This is obviously a regime in which the 
NJL model can not be applied. Our assumption on the coincidence of deconfinement and the chiral phase transitions 
at finite density, has not yet a QCD-motivated argument. Nevertheless, this coincidence has been found in Lattice 
QCD calculations at finite temperature (see Ref. jS^I and references therein) and it has been also adopted in other 
models for the EoS at finite chemical potential as the N JL-inspired model proposed in Ref. (3^ . Interestingly, within 
the Dyson-Schwinger approach, it is possible to define an order parameter for deconfinement at finite density and it 
turns out that the two phase transitions occur simultaneously [5i|. We will examine here both choices for the bag 
constant Bq and B*. 

We remark that another possible scenario has been proposed for the finite density phase transition in Refs. (5^.[53|: 
there is no deconfinement at all at large density, only the chiral phase transition occurs and the quarks are still 
confined. 



A. Results 



In order to compute the EoS needed for compact stars, the pressure, Eq. ([7]), must be minimised with respect to 
the chiral and color superconducting order parameters, ctq and Ac, and therefore six gap equations are obtained: 

1^ = 0, (8) 
. (9) 
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Moreover, local electric and color charge neutrality are met if three other three equations are satisfied, 



ng ^ 1^=0. (12) 

These conditions fix the values of the three corresponding chemical potentials, /iq, /is and /ig- One more equation 
is imposed to compute the vector current expectation value: 



(13) 
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FIG. 3: Same as in Fig.[2]for Gd = 1.2Gs- When a large diquark coupling Gd is considered there are for both values of B two 
phase transitions, first from hadronic matter to 2SC quark matter and then from 2SC quark matter to CFL quark matter. 

After all these quantities are fixed, one arrives at the pressure of quark matter as a function of the quark chemical 
potential only and one can easily build the Maxwell construction by solving the equation: 

PHAiifJ.) = PQAiif^) ~ B. (14) 

As discussed before we use for B two values, Bq and B^. In the upper left panel of Fig. [2] we show the pressure as 
function of the quark chemical potential for hadronic matter (thick line) and quark matter (thin line, the dashed line 
corresponding to the 2SC phase and the continuous line to the CFL phase). Parameters are B = Bq, Gp = Gs and 
Gv — 0. In this case there is a direct transition from hadronic matter to the CFL phase, see also Refs. [l^ll^. The 
corresponding EoS is shown in the right panel of Fig. [21 Notice the large jump of the baryon density of about ~ 0.45 
fm^'' at the onset of the phase transition. We notice that in Ref. [s^l a similar parameter set, B — Bq, Gd = GS 
and Gv = 0, gives a different result: there is just a transition from hadronic matter to the 2SC phase. Apart from 
the different choice for the hadronic matter EoS, the 't Hooft interaction term is neglected in that calculation. We 
obtain the same result, a transition from hadronic matter to the 2SC phase, by choosing if = in our model. The 
't Hooft term in fact, as observed in Refs.fi?!. [53|. pushes the 2SC-CFL phase transition to lower chemical potentials 
rendering the CFL phase the favored quark matter phase also at intermediate densities. 

Let us study how this result changes if we choose B — B^,, which for this set of parameters is i?, = (424:.8MeV)'^, 
slightly smaller than Bq = {A25AMeV)'^. By construction now there is first a transition from the hadronic matter 
to the 2SC matter and then a second transition from 2SC matter to CFL matter, see the corresponding plots in the 
lower panel of Fig. 

We repeated the previous calculation for Gd — 1.2Gs- For this value of the diquark coupling the superconducting 
gap within the CFL phase at ^ — 500 MeV is Aqfl ^ 160 MeV. Notice that our present knowledge of the CFL gap 
concerns only its order of magnitude i.e. ^ 100 MeV; therefore, a CFL gap larger than 100 MeV is not excluded 
and has been considered also in previous papers [1, [s^, [E^ ■ The effect of increasing the diquark coupling on the 
EoS is to decrease the onset of chiral symmetry restoration, from /ig = 358 MeV for G^i = Gs to = 344 MeV 
for Gd = l-2Gs and the 2SC-CFL phase transition onset, from fig = 415 MeV for Gd — Gs to /ig = 386 MeV for 
Gd — l-2Gs- As shown in Fig. [3] for both choices of B there is a double phase transition with increasing baryon 
density [6^ . 

We include now in our calculation the vector meson term and set Gv = O.2G5 (we consider now only Gd = ^-^Gs)- 
The physical effect of this term is a repulsive interaction between quarks which renders the quark matter EoS stiffer. 
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FIG. 4: Same as in Fig. |3] with the inclusion of a repulsive vector term with Gv = 0.2Gs. In the case B = Bo a, direct transition 
from hadronic matter to CFL matter is found. In the case B = Bt the EoS exhibits instead two phase transitions. 

For the case B — Bq, there is a transition from hadronic matter to CFL matter and for the case B — B^, a. double 
phase transition is present, see Fig. [4] 

IV. MASS-RADIUS RELATIONS 

Let us now discuss the corresponding mass-radius relations for compact stars. In Fig. [5] we show the mass-radius 
diagram for the EoSs with Gd = Gs for B — Bq (EoSl) and B ^ B^, (EoS2) indicated by the dotted lines and the 
EoSs with Gd = l-2Gs for B = Bo (EoS3) and B = B^ (EoS4) indicated by the solid hncs. In the first case, EoSl, 
where a transition from hadronic matter to CFL is found (see upper panel of Fig. [2]), the CFL core is unstable and 
therefore the conclusion is that quark matter does not occur at all in compact stars in agreement with the findings of 
Ref. [29] . This is due to the fact that the transition to CFL matter occurs at a large pressure and chemical potential 
which, as demonstrated within the toy model of Sec. I, strongly disfavours stable configurations. Decreasing the 
value of B to B^, which in the toy model would correspond to a change of the transition pressure, changes the result 
significantly: a stable core of quark matter can be present but only in the 2SC phase, the subsequent transition to 
the CFL phase renders the star unstable. This scenario agrees with the one proposed for instance in Ref.js^l- Let us 
discuss now the cases in which Gd = l-2G's- Obviously the larger the diquark coupling the more favoured is quark 
matter with respect to hadronic matter because both the onset of the chiral symmetry restoration and the one for the 
2SC-CFL transition are shifted to lower densities. This implies that the pressure of the phase transition is smaller 
than in the previous cases but also that the jump in density due to the Maxwell construction is larger (65| . Among 
these two effects, the first favors the stability of the star while the second disfavors the stability of the star (see our 
toy-model discussion). We find that the first effect dominates and hybrid stars containing both the 2SC phase and 
the CFL phase are stable in both cases, B — Bq and B — B^, [66J. The structure of these stellar objects is extremely 
interesting: it contains a crust of hadronic matter matter, a layer of 2SC phase and a core of CFL phase. Possible 
astrophysical implications of these stellar compositions have been already discussed in Refs. [2^, |2^, |23] in connection 
with GRBs. The formation of first the 2SC phase and then of the CFL phase, during the evolution of the star, would 
produce two separate neutrino emission and possibly also two different GRBs emission periods for which there are 
already observational hints |26| . 

We point out that the solutions of the TOV equations obtained without vector interactions reach a maximum mass 
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FIG. 5: Mass-radius relations for the EoSs with Gd = Gs (dashed hnes) for B — Bo and B = B^ and with Gd ~ 1.2Gs (sohd 
hnes) for B = Bq and B = B^,. Gv = in aU cases. The thick curve stands for the hadronic stars. In the case Gd = Gs no 
stable CFL cores are found, only a 2SC core is possible for B = B^. In the case Gd = 1.2Gs, hybrid stars having a CFL core 
are stable. Furthermore, in both cases, B — Bo and B — B, a. layer of 2SC phase is present inside the hybrid stars. 
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FIG. 6: Mass-radius relations for the EoSs with Go = 1.2Gs and Gv = 0.2Gs for B ^ Bo and B = B^ (solid lines). The 
thick curve stands for the hadronic stars. In both cases stars having a CFL core are stable (see the magnification for the case 
B = Bo), in the case B — B, a. double phase transition is present. 

of less than IAMq and therefore are excluded by pulsar mass data. For the case Gv — 0.2Gs we obtain hybrid 
stars with a stable core of CFL phase in the case B = Bq and hybrid stars with both the 2SC and the CFL phases 
for _B = as in the previous cases. These solutions have a corresponding maximum mass of I.SMq (see Fig. |6]) 
due to the effect of the repulsion given by the vector interaction and are not ruled out by the presently available 
astrophysical data. For even larger values of Gy, we arrive at stable hybrid stars with just a CFL core. 

A final important remark concerns the large value of the diquark coupling that we must use to obtain stable quark 
matter cores. The crucial quantity for the stability is a low value of the pressure of the phase transition to quark 
matter as seen already on our discussion of toy-model. A low transition pressure can be obtained for a large diquark 
coupling but also for a reduced constituent mass of the strange quark. Within the NJL model discussed here, nig 
turns out to be quite large. Other quark matter studies within the Dyson-Schwinger approach point towards lower 
values of nis and consequently the CFL phase dominates in an enlarged density region [57( . Unfortunately, the EoS of 
quark matter within the Dyson-Schwinger approach is not yet available. Our analysis suggests that the corresponding 
results for the TOV solutions can be very similar to the ones obtained here. 
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V. CONCLUSIONS 



We have studied the phase transition from hadronic matter to quark matter and the corresponding impact on the 
mass-radius relation of compact stars. By using first a toy-model for the quark matter equation of state wc have 
analysed how the stability of hybrid stars depends on the properties of the phase transition. In particular a low value 
of the pressure at the onset of the phase transition seems to be the most crucial quantity for the stability of the 
hybrid star. We have then computed the equation of state of quark matter within the NJL model by including effects 
from the chiral condensates, the diquark coupling pattern and a repulsion vector term. For large enough values of the 
diquark coupling strength, hybrid stars containing a CFL core arc found to be stable. Moreover, exciting stable hybrid 
stars containing both a layer of a 2SC phase and a core of a CFL phase appears to be possible. This result opens 
the possibility to find new signatures of the presence of quark matter in compact stars. The formation of different 
quark phases in compact stars can release a huge amoimt of energy powering both energetic bursts of neutrino and 
gamma-ray signals. A detailed study of the time evolution of a proto-neutron star including the possible appearance 
of the 2SC phase followed by the formation of the CFL phase could reveal some new features and opportunities for 
detecting the chiral phase transition of QCD in the sky. 
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